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Abstract. We present a detailed analysis of the carbon and nitrogen abundances of 
two dwarf carbon-enhanced metal-poor (CEMP) stars: SDSS J1349-0229 and SDSS 
J0912+0216. We also report the oxygen abundance of SDSS J1349-0229. These stars are 
metal-poor, with [Fe/H] < -2.5, and were selected from our ongoing survey of extremely 
metal-poor dwarf candidates from the Sloan Digital Sky Survey (SDSS). The carbon, nitro- 
gen and oxygen abundances rely on molecular lines which form in the outer layers of the 
stellar atmosphere. It is known that convection in metal-poor stars induces very low temper- 
atures which are not predicted by 'classical' ID stellar atmospheres. To obtain the correct 
temperature structure, one needs full 3D hydrodynamical models. Using C05BOLD 3D 
hydrodynamical model atmospheres and the Linfor3D line formation code, molecular lines 
of CH, NH, OH and C2 were computed, and 3D carbon, nitrogen and oxygen abundances 
were determined. The resulting carbon abundances were compared to abundances derived 
using atomic C 1 lines in ID LTE and NLTE. For one star, SDSS J1349-0229, we were able 
to compare the 3D oxygen abundance from OH lines to O 1 lines in ID LTE and NLTE. 
There is not a good agreement between the carbon abundances determined from C2 bands 
and from the CH band, and molecular lines do not agree with the atomic C 1 lines. Although 
this may be partly due to uncertainties in the transition probabilities of the molecular bands 
it certainly has to do with the temperature structure of the outer layers of the adopted model 
atmosphere. In fact the discrepancy between C2 and CH is in opposite directions when us- 
ing 3D and ID models. Confronted with this inconsistency, we explore the influence of the 
3D model properties on the molecular abundance determination. In particular, the choice of 
the number of opacity bins used in the model calculations and its subsequent effects on the 
temperature structure and molecular line formation is discussed. 
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1. Introduction 

Send offprint requests to: N. T. Behara SDSS J1349-0229 and SDSS J0912+0216 

were selected from our ongoing survey of ex- 
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tremely metal-poor star candidates from the 
Sloan Digital Sky Survey. SDSS J 1349-0229 
has an effective temperature of 6200 K, log g 
= 4.00 and [Fe/H] = -3.0. SDSS J0912+0216 
is slightly hotter with a temperature of 6500 K, 
log g = 4.50 and [Fe/H] = -2.5. These two stars 
belong to the class CEMP-r+s, since their at- 
mospheres are enhanced in carbon, and show 
overabundances of s and r-process elements. 
For details on the abundances of the neutron- 
capture elements, see Behara et al. (2009). 

The focus of this study is on the abun- 
dances of CEMP stars determined predomi- 
nately using molecular lines: carbon, nitrogen 
and oxygen. It is known that overcooling in the 
outermost layers of metal-poor stars produced 
by convection is not predicted by classical ID 
stellar atmospheres and is best reproduced by 
3D hydrodynamical simulations (Asplund et 
al. 1999; Caffau & Ludwig 2007; Gonzalez 
Hernandez et al. 2008). An accurate descrip- 
tion of these outer layers is crucial for the de- 
termination of the molecular abundances. 

Fortunately these two CEMP stars are hot 
enough to display atomic carbon lines. These 
lines are formed much deeper in the atmo- 
sphere compared to molecular lines, and are 
thus shielded from the overcooled region of the 
atmosphere, allowing for a validity check of 
the abundances obtained from molecular lines. 

2. Abundance analysis 

In our analysis, we used 3D model atmo- 
spheres computed with the CO 5 BOLD code 
(Freytag et al. 2002; Wedemeyer et al. 2004). 
The 3D spectral synthesis calculations were 
performed with the code Linfor3D. The pa- 
rameters of the 3D models used are listed in 
Table □ 

We compared each of our 3D models to a 
corresponding standard hydrostatic ID model, 
and define the 3D correction of an abundance 
measure in the sense 3D - ID. The ID mod- 
els are calculated using plane-parallel geome- 
try and employ the same equation-of-state and 
opacities as the CO 5 BOLD models. 

With the exception of C2, all of the 
abundances were determined using equivalent 
widths of unblended lines. The C2 abundance 



Table 1. Parameters of the 3D models. 
The model a was used in the analysis of 
SDSS J 1349-0229. We interpolated between 
the corrections from the models b for SDSS 
J0912+0216. The first five models were com- 
puted using 6 opacity bins, while the last three 
using 12 bins - see discussion for details. 



Model 


T c ff 


logg 


[Fe/H] 


d3t63g40mm30n01 fl 


6270 


4.00 


-3.0 


d3t65g45mm20n01* 


6530 


4.50 


-2.0 


d3t65g45mm30n01 6 


6550 


4.50 


-3.0 


d3t63g40mm20n01 


6280 


4.00 


-2.0 


d3t63g40mml0n01 


6260 


4.00 


-1.0 


d3t63g40mm30n02 


6240 


4.00 


-3.0 


d3t63g40mm20n02 


6250 


4.00 


-2.0 


d3t63g40mml0n02 


6250 


4.00 


-1.0 



Table 2. 3D and ID abundances ex- 
pressed with respect to solar abundances of 
log (e) = 8.50, 7.86, and 8.76 for carbon, nitro- 
gen and oxygen, respectively. 



SDSS J1349-0229 


Element 


[X/Fe] 1D 


[X/Fe] 3D 


CH 


2.82 


2.09 


c 2 


3.16 


1.72 


C 1 NLTE 


2.42 


2.51 


NH 


1.60 


0.67 


OH 


1.88 


1.70 


O 1 NLTE 


1.63 


1.69 


SDSS J0912+0216 


CH 


2.17 


1.67 


C 1 NLTE 


1.38 


1.44 


NH 


1.75 


1.07 



was determined by line fitti n g usin g the same 
code used by ICaffau et al.1 [20051) . For CH, 
we adopted the line list used by Bonifacio et 
al. (1998), for the UV OH lines the gf values 
have been computed from the lifetimes calcu- 
lated by Goldman & Gillis (1981), and for C 2 
and NH we used the molecular line lists by 
Kurucz (2005). The results for the two stars are 
listed in Table [2] 



Behara: 3D molecular lines in CEMP stars 



737 




-8 -6 -4 -2 2 4 

loq T m „ 



Fig. 1. The temperature structures of the 3D 6 bin 
model d3t63g40mm30n01 (top) and 12 bin model 
d3t63g40mm30n02 (bottom) are plotted as a func- 
tion of optical depth along with the average 3D 
(solid line) and corresponding ID (dashed line) tem- 
perature structures. Overplotted on the figure are the 
ranges of depth of formation of the C2, CH, C 1, NH, 
OH and O 1 spectral lines used in the analysis. 

3. Discussion 

The sensitivity of the molecular lines to tem- 
perature is clear from the significant 3D correc- 
tions. In Figure Q] we have plotted the temper- 
ature distributions for the models used in the 
analysis of SDSS J1349-0229. The ranges of 
depth of formation of the lines used in the anal- 
ysis are overplotted for comparison. The differ- 
ent depth of formation of the C2 and CH lines 
in the ID and 3D models explains why in ID 
we derive a larger C abundance from the C2 
lines than from the CH lines, while in 3D the 
reverse is true. The C2 lines are formed much 
higher up in the atmosphere compared to the 
CH lines, and this tendency is much larger in 



Table 3. Comparison between 3D abundances 
derived using models with 6 and 12 opacity 
bins. 



Element 


6 bin [X/Fe] 3D 


12 bin [X/Fe] 3D 


CH 


2.09 


2.22 


c 2 


1.72 


1.99 


C 1 NLTE 


2.51 


2.50 


NH 


0.67 


0.87 


OH 


1.69 


1.69 


O 1 NLTE 


1.69 


1.68 



the 3D than in the ID model. The 3D mod- 
els on the other hand do not achieve a better 
consistency between CH and C2 lines. The C 1 
lines are quite insensitive to 3D effects. Neither 
in 3D nor in ID we achieve consistent results 
between molecular lines and C 1 lines. 

To shed light on the discrepancy between 
the carbon abundance indicators, we use our 
best indicator, C 1, as a reference and ex- 
plore the influence of the 3D model proper- 
ties on the molecular abundance determination. 
In particular, we explore the effect of increas- 
ing the number of opacity bins in the opac- 
ity grougsjrom 6 bins to 12 bins as described 
in iLudwig et alJ |2009). The new abundance 
obtained for SDSS J 1349-0229 are listed in 
Tabled 

The temperature distribution of the 12 bin 
model is shown in Figure Q] below the 6 bin 
model. The different binning scheme results 
in a decrease of the overcooling, thus lower- 
ing the 3D corrections of the molecules formed 
higher in the atmosphere. 

The atomic lines, formed deeper in the at- 
mosphere are unaffected. The 12 bin models 
give a better agreement between the CH the 
C2 lines and both molecular features appear to 
be in better agreement with the C 1. Although 
these results appear encouraging, the forma- 
tion region is so external and at low density, 
that the validity of LTE for molecule formation 
and levels population should be questioned. 

To explore the sensitivity of the models to 
the different binning schemes, we compute 3D 
corrections for CH, C 1, NH, OH and O 1 for 
[Fe/H] = -3.0, -2.0 and -1.0 using the models 
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Fig. 2. 3D corrections as a function of metallicity for CH, OH, NH, O i and C i. The 6 bin models 
are plotted in black, while the 12 bin models are plotted in grey. The OH correction have been 
calculated assuming a typical CEMP composition: an enhanced carbon-to-oxygen ratio of +1.0. 
Overplotted as dashed lines are the corrections using scaled solar abundances. 



listed in Table Q] The corrections are shown in 
Figure[2] The OH lines computed using scaled- 
solar abundances appear the most sensitive. In 
carbon-enhanced atmospheres however, the ef- 
fect is much reduced. Overall, with regards to 
the molecular lines, the sensitivity to the bin- 
ning scheme as well as the 3D corrections de- 
crease with increasing metallicity. The atomic 
lines proved to be immune to both the 3D 
model properties and 3D effects. 
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